Wave Packets can Factorize Numbers 
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We draw attention to various aspects of number theory emerging in the time evolution of ele- 
mentary quantum systems with quadratic phases. Such model systems can be realized in actual 
experiments. Our analysis paves the way to a new, promising and effective method to factorize 
numbers. 
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I. INTRODUCTION 

Interference of waves instead of Newton's light cor- 
puscles — what a revolutionary idea of Thomas Young 
pronounced in his famous talk at the Royal Society in 
1801 [01 . The opposite idea — particles as waves — 
proposed by Louis de Brogiie and cast into the appro- 
priate mathematical form by Erwin Schrodinger more 
than 100 years later opened a new realm in atomic 
physics with more applications today than ever. "Can 
a quantum-mechanical description of physical reality 
be considered complete?" This famous question raised 
by Albert Einstein, Boris Podolsky and Nathan Rosen 
in their seminal paper |^ brought out most clearly the 
importance of entanglement |3| as the basic ingredient 
of quantum theory. Interference, wave mechanics and 
entanglement are the three ideas on which Peter Shor 
could build his factorization algorithm Q with a poly- 
nomial rather than exponential speed. This discovery 
gave birth to the new field of quantum information 
Guided by the same three central ideas we show 
in the present paper that the time evolution of elemen- 
tary quantum systems carries intrinsically the tools to 
factorize numbers @. 

Stationary states have played an important role in 
the development of quantum mechanics. For this rea- 
son the time evolution of wave packets has for many 
years taken a back seat even in textbooks on quan- 
tum theory. In this context it is interesting to recall 
that already in 1927 Earle H. Kcnnard Q studied the 
time evolution of wave packets in simple potentials 
Recently, sophisticated techniques to create short 
laser pulses have opened a new era for a controlled 
creation and observation of wave packets in atoms Q , 
molecules ||l^ and cavities [|ll| . Even wave packet dy- 
namics of cold atoms [0 and Bose-Einstein conden- 

These examples illus- 



sates 1 13 has been observed, 
trate convincingly that today's technology allows to 
promote this gedankenexperiment of factorization to 
a real implementation in a quantum system. 

The use of wave packet dynamics to factorize num- 
bers has been recognized for special examples. Ref- 
erence uses a A^-slit Young interference setup. 
Hence, the number to be factorized is encoded in the 
number of slits. The indicator of a factor is the in- 
terference pattern on a screen as a function of the 



wavelength of the light or matter wave. Whenever the 
intensity across the screen is constant, the appropri- 
ately scaled wavelength provides a factor. Similarly, 
Ref. ||l5| considers the time evolution of a particle in 
a box. In contrast to these contributions we provide a 
general framework for factorization and moreover, ex- 
haust the potential of interference by capitalizing on 
properties of quadratic phase factors. 

Our paper is organized as follows. In Sec. II we 
briefly review the essential ingredients of wave packet 
dynamics p^ . In particular, we present an approx- 
imation for the autocorrelation function that illus- 
trates its ability to factorize numbers. The richness of 
the autocorrelation function originates from quadratic 
phases which are also the origin of the Talbot effect 
pj^ , revivals and fractional revivals Jl^ , curlicues 
and quantum carpets pO| . In Sec. p]| we turn to one 
of the most elementary quantum systems, namely the 
quantum rotor. In this case the autocorrelation func- 
tion is a GauB sum which has been discussed exten- 
sively in number theory |^l[. We show that under 
appropriate conditions factors of a number manifest 
themselves in non-vanishing real or imaginary parts 
of the autocorrelation function. We explain this phe- 
nomenon using the properties of Gauf5 sums. Our con- 
clusions and outlook presented in Sec. IV summarize 



the main results and provide an "idea for an idea" how 
to incorporate entanglement into our scheme. 



II. WAVE PACKETS 

In the present section we briefly review the essential 
features of wave packet dynamics. We start from the 
definition of the autocorrelation function and approx- 
imate the energy spectrum by one that is quadratic in 
the quantum number. Physical systems such as Ryd- 
berg atoms and diatomic molecules are well described 
by such a spectrum. We then show the possibility of 
factorizing a number. 



A. Autocorrelation function 

In standard wave packet experiments many phys- 
ical observables of an atom or molecule are accessi- 
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ble. However, one quantity is most appropriate to ob- 
tain information about the dynamics and is nowadays 
measured routinely: It is the so-caUed autocorrelation 
function 

oo 
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which indicates the overlap between the time-evolved 
wave packet 
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and the initial state. Here, the coefhcients ipn denote 
the probability amplitudes of the initial state in the 
basis of energy eigenstates |n) with corresponding en- 
ergy eigenvalues i?„. 

While so far the autocorrelation function has served 
as a tool for observing the dynamics, our goal is to use 
its intrinsic features to perform computational tasks. 
In particular, we show that autocorrelation functions 
can reveal prime factors of an appropriately encoded 
number. 

For this purpose, we consider a wave packet which 
is localized in energy space around a quantum number 
n, that is, the coefhcients ■(/'« assume only significant 
values in a small vicinity of n = fi. This feature allows 
us to expand the energy spectrum 
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around n. Here we have assumed that the spectrum is 
smooth and dense. Moreover, we have omitted higher 
order terms. 

The approximation, Eq. (^, is valid for smooth po- 
tentials and a narrow energy distribution of the wave 
packet. To be specific, we consider the sodium dimer 
Na2. Here a proper laser pulse lifts the vibrational 
ground state in the potential of the electronic state 
X^'^Eg to the excited electronic state A^E^. In the 
potential curve of the excited electronic state the wave 
packet is no longer the ground state, but has a small 
spectral width around a mean value En- Indeed, in 
this case third order terms can be neglected 

To study the wave packet dynamics, Eq. (0), we use 
the approximation, Eq. (|^), and obtain the phase 

Ent/h^ Ent/h+^{n-n) + ^{n-ny. (4) 

Here we have introduced the classical period T^i and 
revival time T. 

Thus the autocorrelation function, Eq. (^, reads 



S{t) 



-iEat/h 



m— — OO 



11 



13 



17 



19 




11 



13 



17 



19 



-0.2 



-0.1 





At 



0.1 



0.2 



FIG. 1: Factorization based on the autocorrelation func- 
tion of a wave packet: a "guess-the-factors riddle". We 
show the absolute value squared |5'jv(r = £ + At)P of the 
autocorrelation function, defined in Eq. (^, as a function 
of dimensionless time At in the vicinity of various integers 
£ — 2,3,5 . . .. The autocorrelation function has a symmet- 
ric maximum at an integer, that is at the origin of each 
horizontal axis, provided this integer is a factor of N. In 
the present case we clearly recognize 5 and 11 as factors of 
A'^ — 55. The width of the initial wave packet is Am = 10. 



(5) 



with a shifted summation index m = n — n and the 
weight function Wm = iV'm+fiP- Moreover, we have 
extended the lower bound of the sum to —oo, since in 
this domain the weight factors Wm do not contribute 
significantly. 
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B. An allusion to factorization 

In order to demonstrate the ability of the autocor- 
relation function to factorize a number, we start from 
a Gaussian weight function Wm oc exp [— m^/2Am^] 
and fix the ratio N = T/T^i to an integer. In this case, 
we have to deal with the sum 



SNir) = ^ WVexp 
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for dimensionless time t = t/Td- 

The sum Eq. is already the key to factorization, 
although it might not be obvious at this moment. In 
order to bring this latent capability to light, we fo- 
cus in Fig. |l| on the modulus squared |S'Ar(r)|^ of the 
autocorrelation function for = 55. We show this 
quantity for times r in the vicinity of prime numbers 
£. Only if £ is a factor of N the function displays a 
clear symmetric maximum. 

The phenomenon of fractional revivals lies at 
the heart of this feature. It is well-known that wave 
packets can show revivals at multiples of the time T 
and fractional revivals at times t = with p and q 
being mutually prime integers. In our case, this fact 
implies that at integer multiples of Td — or in other 
words for integers t — we find a fractional revival 
of the wave packet, that is, a symmetric sequence of 
Gaussian shaped humps. A detailed analysis of this 
structure Q brings to light that in most cases neigh- 
boring Gaussians overlap and interfere. For this rea- 
son, at most times no distinguished revivals appear. 
Only for times r equal to a factor of N the peaks are 
sharp enough so that the central peak does not overlap 
with its neighbors and reveals a clear maximum. 



III. QUANTUM ROTOR 

The crucial point in our factorization scheme is the 
interference of quadratic phase factors. In the previous 
section, they were the consequence of an anharmonic 
potential. But quadratic phase factors even appear 
for a free particle with a continuous energy spectrum. 
However, in order to encode integers we need a dis- 
crete one. Such a spectrum can be achieved, for ex- 
ample, by periodic boundary conditions. A diatomic 
molecule rotating perpendicular to its molecular axis 
is one example of such a quantum rotor. 



A. Time evolution 

In order to lay the foundations for the factoriza- 
tion scheme presented in this section we first briefly 
review the essential properties of the one-dimensional 
quantum rotor. Then we derive the corresponding au- 
tocorrelation function. 



The time evolution of the quantum rotor is governed 
by the Hamiltonian 



H= — 
21 
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with the operator J of the angular momentum and the 
moment of inertia /. Periodic boundary conditions 
enforce a discrete spectrum of J with eigenstates |fc) 
and eigenvalues hk with fc e Z. 

The resulting discrete energy spectrum E^. = 
H'^k'^/(2I) gives rise to the time evolution 
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of the quantum rotor. 

We model an initially well-localized wave packet us- 
ing the expansion coefficients 



, ^ f 1/VN for < A: < iV 
\ otherwise 

which leads to the explicit expression 
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for the autocorrelation function. Here we have intro- 
duced the revival time T = Airl/h. 

Again quadratic phases appear. However, the linear 
term present in the expression Eq. (H) is absent. 



B. Factorization 



In order to bring out most clearly the number N to 
be factorized, we cast the sum, Eq. (10), in the form 
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with n = e Z. We emphasize that only for in- 
tegers n the autocorrelation function Sn is periodic 
with period N . 

In Fig. 11(a) and (d) we show |S'jv(?i)P for n = 1 
to N, with = 21 and N = 15, respectively. We 
recognize dominant maxima for values of n whenever 
n and N share a common factor. The appearance of 
these maxima already indicates that a sum of only 
quadratic phases allows us to find prime factors of N . 

We can even go a step further and concentrate 
on the real and imaginary part of Sn, depicted in 
Fig. 1(b) and (c) for TV = 21, and (e) and (f) for 

= 15, respectively. It is remarkable that for = 21 
only numbers which share a common factor with A^ ap- 
pear in the imaginary part. Fig. |^(c), whereas all the 
other ones vanish. This is not the case for A^ = 15, 
where only one factor and its multiples appear in the 
real part. Fig. |2|(e). 
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FIG. 2: Autocorrelation functions of the quantum rotor for A*' = 21 (left column) and A = 15 (right column). From top 
to bottom, we show modulus squared, real and imaginary part of Sjv, Eq. ([llj). For N = 21 only multiples of both factors 
3 and 7 show non-zero imaginary parts. For A = 15 only multiples of the factor 3 yield a non-zero real part. 



C. Gaufi sums 

We recognize that the autocorrelation function, 
Eq. (O), is of the form of a Gaufi sum 
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m=0 
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Therefore, the structures found in Fig. || result from 
the intrinsic properties of Gaufi sums. In order to 
highlight this connection, we first briefly review the 
properties of Gaufi sums and then apply the results to 
the autocorrelation function Eq. (pd]). 

For odd h and mutually prime a and b the sum can 
be evaluated pl[ and yields 



G(a,6) 



1 /2a 
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Here we have introduced the Jacobi symbol (^) which 
expresses the number-theoretical complexity of the 
quadratic sum and which can take on the values ±1. 

A closer look at Fig. ^ reveals that the emergence 
of factors in either the real or imaginary part depends 
strongly on the representation of the factors as 4s -I- 1 
or 4s -|- 3 with s being a non-negative integer. Indeed, 
all odd numbers r are either of the form 4s -I- 1 or 



4s -|- 3, and therefore, they are members either of the 
set 



Ml = {r\r = As + 1} 



7W3 = {r I r = 4s + 3} . 



(14) 



(15) 



For the case 6 G Ali the sum G(a, b) is purely real. 
Indeed, for 4s + 1 we find 



G cx exp 



I— 4s 
4 



exp[i7rs] = (-1)". (16) 



When b G A^s, that is 6 = 4s + 3, we arrive at 

V 



G cx exp i-(4s-|-2) 



exp 



(17) 



which is purely imaginary. 

Moreover, Eq. (O) implies that the modulus 



1 



|G(a,6)| = ^ 



(18) 



is independent of the value of a, provided that a and 
6 are mutually prime. 
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D. An explanation of the patterns 



factors of TV. 



With the help of these properties of the GauB sum 
we can now relate the structures in the autocorrelation 
function 



SN{n) ^ G{n,N) 



(19) 



apparent in Fig. ^ to the factors p and q of N. Here 
we assume that p and q are the only factors of the 
number N, that is N — p ■ q. Moreover, we stipulate 
that p and q are both odd. For the case of even N one 
can easily extract factors 2. 

According to Eq. ( |l^ ) the modulus squared \S]\[\'^ is 
given by 1/N if n and N do not share a common fac- 



tor. Otherwise, the common factor cancels in Eq. (12) 
and we have b— l/porb—l/q and thus with the 



help of Eq. (P) we have \Sn\^ = l/p or 



_ n\' = l/q, 

respectively. 

We are now in the position to explain the patterns 
in the real and imaginary parts using the classification 
of the prime factors into the sets Mi and Ms- We 
have to distinguish three cases: (i) p,q € Aii, (ii) 
p,q & Ms and (iii) p € Mi,q £ M3 or vice versa. 
The cases (i) and (ii) imply N £ Mi, whereas (iii) 
results in N £ M3. 

We start our discussion with an example for case 
(ii). For iV = 21 = 4 • 5 + 1 = 3 • 7 both factors 
p = 3 = 4- 0-|-3 and q^7 = 4:-l + 3 are elements of 
the set Ms and N £ Mi. If n and N do not share a 
common factor then b = N and the value of the GauB 
sum is real. Otherwise, we can extract a factor and get 
b = q OT b = p, respectively. Consequently, the value 
of the sum is purely imaginary, since p,q £ Ms- Since 
the Jacobi symbol mentioned in Eq. ( p^ ) can only take 
on ±1 the corresponding value oia~n/pova = n/q 
does not change this behavior. 

The same arguments hold true for the example N = 
15 = 4 • 3 + 3 G Ms illustrating case (iii). Indeed, the 
factorsp = 3 = 4-0-1-3 £ Ms and g = 5 = 4-1-1-1 £ Mi 
are not members of the same set. Therefore, if n is 
a multiple of p we have 6 = 5 and thus the sum is 
purely real. In all other cases, we have either b = N 
or b = p which are both members of A^3 and yield 
purely imaginary results. 

The last case (i), not shown in Fig. |^, discusses fac- 
tors p, q and N £ Ml, resulting in purely real values 
of the sum for all n. Unfortunately, in this case we 
do not have vanishing values for integers that are not 



IV. CONCLUSIONS AND OUTLOOK 

The time evolution of a wave packet brings to light 
key aspects of number theory. This feature allows us 
to find prime factors p and q of an appropriately en- 
coded number N — p-q. These properties are intrinsic 
to the autocorrelation function, which can be easily 
measured. At the heart of this phenomenon lie the 
mathematical properties of GauB sums. 

Wave packets of the quantum rotor exhibit these 
number-theoretical properties in the cleanest form. 
The most remarkable result is, that in the case of 
p,q £ Ms we find all factors (and its multiples) of 
N = p ■ q in the imaginary part of the autocorrelation 
function. All other integers yield a vanishing imagi- 
nary part. 

But how to make use of these unique features? 
Imagine that Im(S'Ar) represents the statistics of an 
observable O. Thus, each single measurement of O 
would provide one factor of A'^ in a non-probabilistic 
way. Similarly, in the mixed case p £ Mi and q £ Ms 
the measurement of the real part of Sn would provide 
one factor with high probability. 

These insights pave the way to an effective method 
of factorizing numbers. The next step is to find a com- 
posite system which (i) intrinsically includes quadratic 
phase factors (ii) allows parallel calculation of all pos- 
sible factors and (iii) enables us to find a suitable ob- 
servable, which reveals all the information about the 
factors in a deterministic way by use of entanglement. 
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